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Abstract 

All experiments with ultracold atoms are performed in the presence of background residual gas. 
With the help of a suitable master equation we investigate a role of these fast atoms on the loss 
of coherence in optical lattices. We present an exact solution of the master equation and give the 
analytic formulas for all correlation functions in the presence of one body losses. Additionally we 
discuss existing of a Schrodinger cat state predicted in this system in 
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I. INTRODUCTION 



Ultracold bosonic atoms placed in wells of the optical dipole potential formed by the 
interference of counter-propagating laser beams are studied in many recent experiments 
y, y, y, 01- a Mott insulator state t2|l, a collapse and revivals of the coherence in the 

n n 

superfiuid state pLi], a massive entanglement between atoms j6| and an Anderson localization 
were all achieved. In a recent paper [8] a single site addressability has been demonstrated 
- a property essential for the future application of atoms in optical lattices in the quantum 
information processes. In all these processes a potentially significant role is played by losses. 
As always, the losses generate decoherence, a deterioration of the ability to interfere. In 
particular even the best achievable vacuum in the cell holding the ultracold atoms is never 
perfect. There are background atoms that are very fast. From time to time they collide with 
the relevant bosons that are practically at rest. As a result the atom trapped in the well 
is rapidly ejected and flies away. We are losing one atom but also the phase of the system 
undergoes a jolt. As a result, the atomic interference gets smeared-out. In this paper we are 
studying in some detail this phenomenon. To this end we use a suitable master equation. 
We present an analytical solution for the evolution of a large class of correlation functions. 
Then, we use these solutions to discuss the first order coherence for the model of a perfect 
Mott insulator initial state and a model of a superfiuid initial state. In the latter case we 
show that often used the coherent state as a representation of the state in a single well is 
sick. While correctly representing the distribution of the number of atoms found in a single 
site, it incorrectlybreaks the U{1) symmetry of this state. We point out that discussed by 
some authors [ll,!^ Schrodinger cat state is an artifact of the coherent state representation. 



II. MODEL 



A motivation of our work is an experimental demonstration of a collapse and revivals of 
matter waves described in [Ij. The authors of this paper prepared bosonic ^''Rb atoms in 
the superfiuid state in the optical lattice. Suddenly they increased the intensity of lasers 
creating the deep optical lattice, such that the tunneling between lattice sites was completely 
suppressed. They hold the system during time t and then switch-off the trapping potential, 
allowing cold atoms to expand. Next they take a picture of interfering matter waves. They 
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were repeating the experiment extending only the holding time t. In the first few photos 
they recognized an interference pattern typical for a superfluid state [2]. Then this pattern 
had vanished just to a Gaussian noise, but at the time t = t^ev it appeared again. We study 
the sensitivity of these revivals to decoherence due to non-perfect vacuum. Atoms from the 
residual thermal gas collide with ultracold atoms. We assume that the typical energy of an 
incoming atom from the background is so high that during the collision the condensed atom 
is kicked out from the lattice. This kind of one body losses leads inevitably to a reduced 
coherence in the system. 

In order to describe the evolution of cold atoms in an optical lattice that includes these 
collisions we derived the following master equation: 

K 



dtp 
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k=l 



where annihilates one boson in the kth lattice site, K is the number of lattice sites, 7 is 
a damping coefficient, and 



K 



K 



Ha = f\^^o'^a\ak + —'^ala^i.akak, (2) 

k=l k=l 

is the Bose-Hubbard Hamiltonian for atoms in an optical lattice in the absence of tunneling. 
The parameter Tiloq in ([2]) describes the kinetic and potential energy of one atom in the 
lattice, and U is an average interaction energy per a pair of atoms. The details of the 
Hamiltonian Ha are given in [l^. We estimate the damping coefficient 7 and sketch the 
derivation of ([1]) in the appendix. Having analytical form of all coefficients in this model we 
are able to compare our theoretical results of this papers with the experimental papers. 

We assume that atoms in each well evolve independently so the final master equation is 
a simple extension of the master equation for a single well, examined more generally in 
This type of master equation ([T]) was also a starting point in 12 , isl, 14 1 . 

In this paper we expand all states in the Fock basis, where | ni, ■ ■ ■ , n^), means n\ atoms 
in the first well, 77-2 atoms in the second, etc. Here are three interesting states: 

• The superfluid state, called also multinomial state, with particles in K potential 
wells 



ni,n2, ■ ■ ■ ,nK). 



(3) 



rii,--- ,nx 



The sum in expression (|3]) is over all non-negative integers, which adds-up to the initial 



N\ 



number of atoms A^, and ( ) = 

' Vni,--- ,nR-7'-| n\'.---nK' 

in fact a SU{K) coherent state |l5l ]. 



is the multinomial symbol. This state is 



The Mott insulator state, which describes the system with exactly n atoms in each 
well: 

K 



(4) 



The coherent state 



K 



k=l 



k=l \n=0 V^^' J 

Although the coherent state violates the barionic superselection rule, it seems to work 
quite well as an approximation of the superfluid with the total number of atoms 
= K \c (f- The evolution of this state in a context of optical lattices was considered 
i.e. in [iSj. 



III. GENERATING FUNCTION 

The master equation ([1]) describes a system of K independent modes, with the interaction 
in each mode and with the one body losses. In the next section IIVI we use this model to 
describe the influence of background atoms for coherence properties in optical lattices. We 
stress that equation ([T]) has quite general form and might be appropriate for other systems. 
In this section we present the exact solution of a master equation of this type. 

Let us consider a dynamics of an average value of the normal ordered operator: 

(n(4)'*n4'). (6) 

\fc=i k=i I 

At this point we assume for simplicity that for all indices /c, coefficient or Ik is equal 
to 0. It means + h = max{rk, 1^). 

We focus on the most important case of expression ([6]), correlation functions, where the 
number of creation operators and annihilation operators is the same. It means Xl^i ^fc — 
= R, where R is the order of correlation function. Below we will often use this 

property. 
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The density operator has the following matrix form in the Fock basis: 

N 

P= I "-1, ■ ■ ■ , nK){mi, ■■■ ,771x1 . (7) 

ni,--- ,nx 
mi,--- ,mii 

The density matrix elements which are appropriate to calculate ([6]) are of the form 
Pnl+h'-'-' uk+Ik ^ where are non-negative integer numbers. We introduce a shorter no- 
tation for indices {rii + /i, ■ ■ ■ ,nK + J-k} = n + I. Thus / belongs to and with we 
denote the versor in the kth direction in this space {{£k)i = ^k,i)- From the master equation 
([1]) we know their evolutions 



dtPlX: = (rfc - h) i2n, + r, + - 1) plX\ + 

n+r+e, 



+27 Ef ^{nk + l){nk + ru + lk + l)pT"^"' 



-7Ef(2nfc + r, + 4)p-:r, (8) 

where A = U/{2n). 

We define the generating function /ip 

n k=l m=l 

where x is any vector in space. The function in ([9]) is defined in such a way that 
along a line Xi = ■ ■ ■ = = 1 is just equal to Q). Multiplying the equation ([8]) by 
Y[k=i ^k'' YlmXi' V^A: + ^'^d then taking the sum over all indices we derive an equation 
for hi 

K 



dthl = ^ (27 - (27 - 2iA (rfc - Z^)) x^) ^^../i^ + 

k=l 

K 

+ E (^^ (^^ - - 27^) hi 



k=l 

This is the first order partial differential equation, so one can find the solution by the method 
of characteristics: 

xhl{f{xi,n -h,t),--- ,f{xK,rK-lK,t),0), (10) 
where / (x, r, t) = — — and = 27 — 2iAr. 



From the generating function one can easily extract all density matrix terms: 



n+r 
h'n+l 



Solutions (fTOj) and flTTl) are valic 
exact solutions of papers 



nfc=i \/i.^k + 4-)' {^k + rfc)! 



(11) 
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16, 



or all initial states and any number of wells. They cover 



17| as special cases. The paper |16| treats all correlation 
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17l | are interesting discussions 



functions in optical lattices but without losses. The papers 
of one, two and three body losses in two mode system. There are presented exact solution 
for low order correlation function in the case of one body losses using the stochastic wave 
approximation. 

The equation (fTTj) says that the terms of density matrix decay exponentially with the rate 
growing as we go away from the diagonal. For long times all off-diagonal terms vanish and 
the state of the system becomes completely mixed, in agreement with the general properties 
of the decoherence. 



IV. DECOHERENCE 

In this section we will use the general solution (iTTl) for the three states introduced in the 
section IIIII Subscripts 'coh', 'sf and 'mi' denote the coherent state, the superfluid state 
and the Mott insulator, respectively. For shorter notation we omit indices r and Z, namely 
h\ = h. The appropriate generating functions have forms 

K 



K X N~R 



X 



K 



5^/(xfc,rfc-4,t) (12) 

hcoh{x-t) = JJ («*)'^'= («)''= el"l'(^("'=''''=-''='*)-^) (13) 
fc=i 

K 

h^dx-t) = J](l-(l-a:fce-27*))^/^ (14) 

k=l 

Using these functions we can analyze all correlation functions for these three different initial 
states. It is not surprising that for all three states the average number of atoms per well is 
given by the same formula 

(ala,) = 9,,/^°(l,---,l;t) = ^e-2^*. (15) 



This exponential decay of the particle number was observed in several experiments. The 
coefficient (27)^^ is known in literature as the lifetime of the trap. In the appendix we give 
the straightforward but lengthy derivations of 7: 

7 = cr cross Pb{v) /VQ, 

where (7 cross IS Qb cross section for the collision between 'hot' and cold atoms, ph is a density 
of the background atoms and (v) is their average velocity. We compare measured values 
of 7 with our analytical estimation getting quantitative agreement p^, 1191, 120|, except the 



experiments with the metastable helium *He 



21 



221 . We suppose that in experiments with 



helium the interaction with background atoms is not the main source of one body losses. 



A. First order correlation function 



More interesting is the evolution of the first order correlation function 



9' 



{a\aj) 



{a]ai) {a]a-j) 



n 

which in the experiment of Bloch [1], is responsible for the revivals and collapses of the 
interference pattern. For the Mott insulator state this function is just equal to - to describe 
interference of matter waves in this state higher order correlation function are needed. If 
the initial state is a superfluid or a coherent state one get 

91h = el-l''"^^ (17) 

where w{t) = 2 (^-^{^^-^^^^^'^^^ _ { 
and large number of atoms, both first order correlation functions are identical and equal 
to = 1 + \afw{t). Thus, in typical experiments with optical lattice the damping of 
coherence is the same for superfluid and coherent state. It suggests that the coherent state 
can be used instead of multinomial one also for studying the decoherence in optical lattice, 
but in the next section we are showing the limitation of this replacement. 

We present these functions in figure [H Indeed the evolution of the first order correlation 
function has a structure of collapses and revivals at the multiples of the revival time trev = f ■ 



In the limit of large number of wells 
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FIG. 1: (Color online)The first order correlation function for systems initially in superfluid or 
coherent state. Solid red line, - the damping like in recent experiments 7 = 0.001/s, green dashed 
line - with much stronger damping 7 = 10.0/s. Number of wells 8000 and number of atoms 20000 
- like in the experiment 

2 

For experimental value of parameter 7 no effects of decoherence can be seen. It means that 
in the experiment [l^ there are more important sources of decoherence. As authors of this 
paper suggest it can be the noise of the laser beams. 

We investigate the decoherence varying the experimental parameters - the number of 
atoms, the number of potential wells and the damping parameter 7, but holding constant 
the interaction parameter U. We observe that for bigger ratio K/N revivals are broader and 
collapses of coherence are not complete (see figure [2^). In it case even for large damping, the 
function tends nor to 0, but to some positive value. In the case of strong damping it is 
difficult to determine if the master equation is still valid - it should be verified experimentally. 
On the other hand if we increase the total number of atoms, the revival becomes extremely 
narrow and extremely hard to detect. It is worth to stress again, that we hold the parameter 
U constant and equal to the value from the experiment [1]. In fact U can be changed so 
the time scale might be totally different (like in typical experiments with a double well 
potential). 

We present this two regimes in the figure [2l One can easily see that the first regime, 
where the number of potential wells is bigger than the number of atoms is less susceptible 
to decoherence. Increasing the number of atoms leads to a fast decay of coherence. 
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b)N<K 




FIG. 2: (Color online)Two different regimes for evolution first order correlation function: a) 
number of atoms 1000 and wells: solid red line 1000, dashed green line 10000. b)number of wells 
1000 and atoms: green line 10000, red line 100000. Damping coefficient in all cases 7 = 10/s. 



V. Q-FUNCTIONS DEATH OF A SCHRODINGER CAT 

Up to now we haven't seen any notable consequences of using coherent states for the 
superfluid state. Despite the coherent state in the context of matter wave violates the fun- 
damental law of nature, the barionic number superselection, it seems to be indistinguishable 
from the superfluid state. In the kind of thermodynamic limit the distribution of number of 
atoms per one well is identical, the damping of number of particle is the same, and the first 
and the second order correlation functions are so similar in both cases that in any experi- 
ment they cannot be distinguished [l6j . On the other hand the coherent state is much easier 
to handle than the multinomial one. Hence it is no wonder that instead of multinomial 
state usually only its coherent counterpart is investigated. In its evolution at the time ^trev 
surprisingly appears a Schrodinger catlike state - the state of the system is a superposition 
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of two coherent states with opposite phases . It was even predicted in that this super- 
position state may stable even in the presence of one and three body losses. To shed some 
light on this issue we compute one of quasiproababilty functions. The simplest here is the 
Husimi Q-function: 

Qp (/3) = ^(/3 I p I /3), 

where /3 is a vector in K dimensional complex space and | /3) = (S)a^i I l^k) is a product 
of coherent states in each well. 

We use the generating functions f|T3|) . f|T^ and f|T2l) and the general form for the denisty 
matrix (ITT]) to calculate the Husimi functions for our three initial states in the simplest 
7 = case: 

^ K oo 

QcoH{(3;t) = ^HE ■ ^' X 

k=l n,m 



n\m\ 



Qsf (/3;t) 



K 



N 



N 

ni, ■ ■ ■ ,nfc 

21^/ K , n i2n 



K 



AXtn', 



k=l 



k=l 



{N/K)\ 



The immediate conclusion is that the Q-functions for the coherent and the superfluid 
state are quite different. The first one at time t = is a product of Gaussian functions 
around points a, while the second function depends only on the absolute value of a, so it 
is symmetric around 0. To visualize the Husimi functions, we integrate them over {K — 1) 
lattice sites 

Q{(3;t) = j d/?2---d/3xQ(/5,/?2,---,/3x;t) 

These reduced Husimi functions are presented in figures [3] and HI The reduced Husimi 
function of the superfluid and the Mott insulator state is time independent. 

Thus, we conclude that the much discussed Schrodinger cat for atoms in the optical lattice 
is an artifact of the coherent state representation of the superfluid state. Our conclusion 
would be the same if we used the Wigner function rather than the Husimi function. 
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a) t=0 




FIG. 3: (Color online)The Husimi function Q{t) at times a) 0, b) 2^rev if the initial state is 
coherent. At the time ^trev the system is a superposition of two coherent states. 

VI. CONCLUSIONS 

The most important result of this paper is the exact solution for the master equation ([1]). 
We use this solution to discuss the influence of 'hot' background atoms for the coherence 
properties of ultracold atoms in optical lattices. The calculated damping of the average 
number of particles agrees with experiments. In the main experiment under consideration 
[l| the damping of due to collisions with the background atoms in negligible. We consider 
additional isotropic harmonic confinement also. Then in each potential well atoms have an 
additional energy shift which value depend on the position of the site. The spread of the 
shift's values is however so small, that the effects related to them are negligible. Studying 
further this model we find interesting regimes in which the function does not collapse, 
even for large losses. It means that the interference pattern should be seen in this regime at 
any time so the decoherence is quenched. Furthermore we have shown that using a coherent 
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Superfluid 




Mott Insulator 




FIG. 4: (Color online)The Husimi functions for a) the Mott insulator state b) the superfluid state. 
The Husimi function in both cases is time independent and has U(l) symmetry. 

state for description of cold atoms might lead to wrong conclusions. The Schrodinger cat is 
an example. 
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APPENDIX: COEFFICIENT 7 

We present an estimation of the coefficient 7 in the master equation ([T]). To calculate 
this parameter one must follow whole derivation of the master equation. The origin of 7 lies 
in the second term in the Hamiltonian: 
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H 



d^rij^ (r) HoiJ (^) + ^ / dVdV'V^^ (r) ^j^ (r') V {r - r') (r') ^/^ (r) + ■ ■ ■ , (A.l) 



where an operator ^jJ'^ (r) creates a bosonic field, T-Cq = — -^^2™" ~'~ ^(^) consists of the kinetic 
energy and a periodic potential and V {r — r') is a potential of interaction between particles. 
We represent the field operator ip (r) as a superposition of BEC and background gas modes 

V / 

where a annihilates one boson in a ground state, (r) is a spatial representation of the 
ground state fo Tio, annihilates one background atom with wave-vector k and an- 
nihilates one boson with wave-vector k, which has been kicked out from the condensate. 
All non-condensed bosons are modeled by the plane waves in quantization volume V. We 
introduce two annihilation bosonic operator dk and c^, because in general, ultracold atoms 
and background gas might be formed from different species. 

At this point we make the crudest simplification, we replace V {r — r') with the contact 
pseudo-potential '^^^"' S (r — r'), where a is th s-wave scattering length, and m is a reduced 
mass of colliding particles. This approximation works quite well for collision of two ultracold 
atoms and it is not justified for a collision between one cold and one 'hot' background atom. 
Still, this simplification offers some insight into this process. We expect that the final result 
is at least of the same order of magnitude as the accurate damping constant 7. 

We limit our consideration only to the first two terms in (]A.1|) . Substitution (\A.2\i into 



f lA.ip leads to plenty of terms. Most of them have no influence on the values of 7. The 



significant part is given by the Bose-Hubbard Hamiltonian 

H = Ha + Henv + Hint- 

It consists of three terms. 

• The free evolution of ultracold gas in an optical lattice. 

Ha = fujJod)d + ^d)d)dd, 

where Tiujq is the energy of a single atom in an optical lattice, U is an interaction 
coefficient 
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The term which describes the evolution of free 'hot' atoms He 



^ Y.k (clck + (ifc^fc) ' where hwk 
of mass m and momentum hk. 



2m 



is a kinetic energy of a free particle 



The interaction between cold and 'hot' atoms 
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H 



int 



2^3/2 



fci,fc2,fc3 

lki,k2,k3ciJi^Ck:,a + h.c 



(A.3) 



where g 



. For shorter notation we have introduced in (1A.3P the operator ^ki,k2M- 



We follow the derivation of master equation based on [23|]. Our derivation is simplified in 
comparison with ll|] . We quote it here to obtain the expression for the damping coefficient 
7. We write the von Neumann equation in the interaction picture 



dp 
'dt 



I 

h 



Hint (t) , P , 

where with a tilde we denote the transformation 

Q ^ ^i{HA+HE)t/hQ^-i[HA+HE)t/h 

Using the von Neumann equation one can easily write the first two terms of perturbation 
series of the formal solution of the density matrix (Born approximation): 



dp 



I 

h 



Hint {t) , P (0) 









Hint (t) , 



(A.4) 



In the next steps we trace both sides over the environment and transform it back to the 
Schrodinger picture. We assume that initially the density operator of the system has a form 

P = Pbec ® Pe, 

where pbec is a density operator of BEG and pe is a density operator of the environment. 
In the derivation we also assume that the change of pbec during time t is negligible and 
the evolution is well approximated by the Markov process. All this operation one can find 
in the following formula for coefficient 7: 



japBEca 



j dT%^^k2Mlk{,k'^,k'.,^PBEca^e'^'' *^^Tr |ci,^Cfc3dfc2P£;(iJ., c[,,^Cfc^ | , (A. 5) 



where ^ = U /h a^d + tuo + ojk^ — 00^2 — ^^ka 
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We assume that the background atoms are in a thermal state and the atoms from d-modes 
are initially in the vacuum state |0)(0|. Thus 

Tr [cicuJkjE\c{,cy^^ = (4iCfci) (l + (cigCfcs)) ^tl^tA 

For high temperature the average occupation of any energy level, (c^Cfc), is much less 
than unity what justifies the approximation: 

After that the formula for 7 is reduced to: 



lapBEca^ = -^ll,,k^,k3^pBEca^e'^^ *^^(4iCfci), (A.6) 



The operator 7fci,fc2,fc3 = 1/ '^^'^ Xlfci fc2 fca *^ '"^ '"^''^ contains Fourier transform of 

(r). We approximate this function by a Gaussian with the width a: 



27rcr^ 



3/4 2 

e 4^ 



Then the operator 7 is simplified to 7fci,fc2,fc3 = — 2 — Sfci h^H^ " ''^^ '^^ '"^^ 
We use the next approximations, typical for derivation of a master equation 



/ 



Stt^ 

k 



j dre*"^ = Ti5 (a) 



where is a principal value. After using these three approximations the expression flA.6p 
has a form 

'^apBEco^ = ~ ^^2^27r)^ / ^^^^ I ^^^^ o-Pbeco^t^^ {c^^Ck^), (A.7) 
In (IA.7P we have omitted the term with a principal value. This term should give infinite 



shift to the energy levels, but does not change the damping coefficient. It appears in many 
calculations, for example in description of a spontaneous emission - after the renormalization 
it represent the Lamb shift. We don't want to discuss this term - our aim is only a very 
rough estimation of the damping coefficient. 
In the expression (IA.7P appears the operator 
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A = tuo + Aa'a 1 feifea 

2m m 

We assume that background atoms have much higher energy than cold atoms. Further- 
more in typical experiments with optical lattices in each lattice site the average number of 
atoms is small, so the contribution of interaction energy is negligible. Finally with good 
approximation we have 

^ ~ ^^^3 _|_ Jk}^ = /\ 
2m m 



Then 6 (-^ + ^k,k,) = -^5 [cose - 
where 6 is an angle between the wave- vectors ki and k^. 



After integrating over k^, the equation (1A.7I) has a form 



7 = J / ^kiki{cl Ck,), (A.8) 

4/i3(27r)Vo ' 

We calculated the quantity (c^Cfc) from the grand canonical ensemble using the bound (n) = 
Sfc(^fc^fc) ~ where is the total number of background atoms, and get 



(4cfc) = — 



V ylrnksT , 

Finally we can integrating over wave-vectors of incoming background atoms and simplify 
equation (lA.Sp : 

7 = s^T^a'pM, (A.9) 

where p\, = N/V is a density of background atoms and (v) is an average velocity of them. 
The last formula up to a constant can be deduced classically. 
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